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INTRODUCTION 


The problem of wave diffraction by a conducting half plane is a classical one, and its 
first nomasymptotic solution was given by Sommerfeld in J 890 lib The exact result for the 
surface current induced on t he half plane involvs Fresnel integral functions , which although 
straightforward to confute, are inconvenient for quick calculations, Hence, there have been 
attempts to use approximate expressions for t he surface currents. Since the surface current 
has a singularity at the edge, it is most natural to try to use a line current at the edge as 
an approximation for the diffraction current function (true current minus the geometrical 
optics current). The accuracy is, however, not too good, In fact, since a line current 
radiates 01 uni directionally in the plane transverse to the line, the radiation pattern of 
the diffraction current is only poorly approximated, The common technique in improving 
the accuracy of the GTD is to use nonphysical, so-called equivalent edge currents, which 
depend on the point where the field is to be calculated, These techniques are compared in 
the reference (2) but it is out of scope to delve into any comparisons in this short study. 
However, the equivalent currents can obviously be expressed in terms of a multi filar series- 
expression (soon to be defined), the first term of which is the monofile, or line c jrent, 
The problem studied here is to find an optimal location for the monofile such that the next 
term has the least effect, to make the approximation through the monofile most efficient. 
The idea lias its origin in the corresponding three-dimensional approximation of scatterers 
by multipole series, [3], [4], 


lUULTIFILAR EXPANSION OF SOURCES 


In references |3], [4] the mujtipole expansion in complex space has been applied to 
source approximation in three dimensions. The theory can be modified to two dimensions 
so that instead of point like mult ipoles, threadlike multifilar elements are used. The problem 
need not be truly two dimensional, since the multifilar currents need not be constant. 

Let us consider a current source J(r) and its vector potential A(r) function 


A(r) ~ p J G(D)j(r l )dV'. 


Here, D is the distance function 


(D 


D zz y(r ~ f r ) * (r - r f ) (2) 

and it can be complex if r f is given complex values, 0(D) is the Green function 

r jkD 


0(D) = 


4nD 


(3) 


If the source is almost axial, i,e., of small extension in th direction transverse to the 
~ axis, and we interested in approximations giving good accuracy for the far field 
calculation, the Green function can be approximated in two • dimensions [x ,y) through the 
Taylor expansion: 

O(D) = G(p) -p' ■ VG(p) + \p'p ' : VVG(p) - . 


(4) 


Here, 0(p) is also a function of z s- . although not explicitly shown* When substituted 
in (1 ), after partial integration, we have the expansion 

A{t) = p (/ 0 6V) + h • Wp) + \% ■ VVG(p) + ...) dz'. (5) 

This is equivalent with the following two-dimensional multipole or multifile expansion for 
the current function: 

J(r) = Mz)i{p) + Uz) • V6(p) + ^(s) : VVS(p) + ... (6) 

The first (zeroth order) term corresponds to a monofihr current component, the next 
one (first order term) a bi filar r the following one a guadrifilar component and so on. The 
expansion (6) follows a similar pattern as given in reference [5] for three-dimensional sources 
and it can be formally obtained directly by expanding the function 6(p — p ') within the 
integral J(f) = J(p,z) = f s 6(p - p*)3 (p f } z)dS f . In (5) and (6), the tf . s are polyadic 
2 -dependent current moments defined by 

I 0 = J s j(f)dS, (7) 

h = j $ p'mds, ( 8 ) 

% = J s ppj(f)dS, (9) 

and the series (6) could be further extended [3], [4]. If the axis of the multifile is moved 
from the z axis to an axis going through a point p = a, instead of (6), we can write, 
suppressing the z dependence, 

J(p) = P 0 (a)6(p ~a)+ 7 ,(«) • Vt(p - 6) + ^7(5) : VVS(p - a) + ..., (10) 

with the new current moments 

7 0 (a) = J J(r)dS = 7 0 , (11) 

7,(5.) = j(p - a)J(f)dS = 7,- al 0 , (12) 

72(“) = J (P ~ o-)(f> ~ 5)J(r)dS = ?2 — al\ - ^ 6,ac, 4- Saif), (13) 

In (13), we have applied the notation I\ = F° r sources sufficiently small in the 

transverse plane and of sufficiently slow variation, the series (5) for the radiation field 
converges quickly and the first few terms will give a good approximation. To find a good 
first- term approximation, we can try to move our multifilar current to a point a which 
makes the first term most dominating. This amounts to finding a point a such that the 
first moment in (12) is minimized. Following the method given in [3], [4], assuming Po ^ 0, 
we find 


3 


( 14 ) 


h’ 

h-io 


h 


]u this paper, we study currents that have constant polarization, Thus, the vector 
potential will also be of the same constant polarization, Denoting this polarization by a 

unit vector u\ we can write for 7() 7(jie,/j - J \u\ etc. Thus, the rank of the current 
moment polyadics can be lowered by one, which makes the formulas appear simpler. For, 
example, (14) can be written us 



( 1 4a) 


Also, it is easy to see, that the choice (14a) not only minimizes the norm of (12), but 
actually makes i\(h) — 0, Further, if the zeroth moment is zero, the following value for a 
will render the first moment the most dominant one by making the second moment (13) 
vanish: 


a 


hi] 

\h \ 2 



(15) 


For a current depending on a single coordinate «,*, (15) can be reduced to the simple 
formula 



(15«) ' 


If a planar surface current is replaced by the multi filar expansion, it can be easily seen 
that the monofilar term gives the far field radiation correctly in the direction perpendicular 
to the plane. By choosing the location of the monofile as given in (14a), it can further 
be shown, that also the derivative of the radiation pattern becomes correct in the same 
direction, To see this, we consider the far field in the direction of the unit vector fi r , as 
A(r) - ftG(r)f{u r ) , with the radiation function 


/(fir) “ J J(r')dV'. (1C) 

If h is the unit vector normal to the plane, obviously n * r 0 when r is on the plane, and 
the origin lies on the plane, For a direction u r differing from n only a little, the term in 
the exponent in (16) is small and we may write 

f(v r ) ^ jj 1 4 - jkur ; r f )J(r')dS f 

* h 4 jku t • 7j » (1 4 jkiir * a)J 0 * (17) 

The last expression is equal to the radiation function due to a monofilar current at the 
location p a. Thus, it gives a two-term approximation for the radiation pattern at the 
direction u r n, or both the radiation and its angular derivative correctly* The monofile 
at p sss 0 gives only the first term correctly. 
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THE HALF-PLANS CURRENT PROBLEM 

Lot us consider first the simple case of plane wave coming perpcndi< ularly (0 q “ 0 and 
00 ~ 0...180 0 , Fig.l) to the edge with the field polarized so that either the electric or the 
magnetic field is in the direction of the edge (z direction)* These two cases are sometimes 
called * respectively* the E wave and the H wave, We prefer calling them the TM wave 
and the TE wave, where "F denotes transverse to the direction of the edge (z coordinate)* 
For the induced surface current we can write 

■M*) - Jgo(r) -I ./,/(*), (18) 

where the geometrical optics surface current is (n ~ Uy) 

Jgo(x ) = 2Uy X H iC i^CO,^ 0 ( 19 ) 

According to the reference |6, p*151], the diffraction current can be expressed for the TM 
Wave as 

jJ M (x) = -2 J™(x) ^F-(r) - , (20) 


and for the TE wave as 


3j E (x) = -2j™(x)J iF-(v). 

(21) 

Here, F_(t>) is a Fresnel integral function and with 


if ss ^/kx{l + COS(f> o), 

(22) 

Jo 

(23) 

and with positive branch for the square root in (22). The Fresnel integral has 

the properties 


(24) 

Flu) ^ (l 1 "V (2ro ” 1 ^ m V 
^ 2 jv y ““(w - l)!(2v ) 2 ' ? * ) ‘ 

(25) 


(25) converges quickly for large v values. It is easy to see that* for large v, the function 
(25), and hence (21), decay as 1/v/?, whereas the function in brackets in (20) obeys the 
law l/xy/x . Hence, the diffraction current in the TM case is more concentrated near the 



LINK (TRRKNT APPROXIMATIONS 


For tilt 1 approximation of the previous surface currents (20), (21), we have to eval- 
nate the* current moments (7), (HJ, which for the TM wave are of the form /q 


n /7*-V P ^ tl M J r . „ a 1 ritw .......... ..r il.~ 'f 7 h 

" - Ml * 7 

; > The moments can be evaluated Invoking; t he following integral expressionsj 


and for the TK wave of the form /o 
uated in 

fOL t V* 


urir-j™ 


L 


0 V* 


dp 


P ?" 4 P*dr * 
Jo 


j; 


• 7T 

V v' 

v/?r( 2n - J )! 

(n 1)!2 2 " 1 p’ H 2 ’ 

-SftvS. ,4). 


( 20 ) 

(27) 


Jo ” 2 

'I Inis, we can write for the mouofilar approximation of the surface diffraction currents 


(28) 


iPw * - « 7 ">. 


«/J /v ( vr ) ss Auy y Ifjljj tl b(p - a 1 1: ). 


7/^ 


„7£\ 


( 20 ) 

(30) 


After a number of algebraic operations, the normalized current moments can be found 
in the following forms 


/TM 

J 0 


yVr xvn&Q 

2jWty \/i f erw^of 3 4 \/3 4 co. 


*0 


1 - ^ 

2A 3 * v'i * eofl0o 


rVW 


v'?f smon 

ly/A* v 1 * msfjntl \ I cosC*(ti" 


r p: 


1 2 4 yt '4 CO50>0 

4P (I t 'vl - ^o) 2 ' 


(31) 

(32) 
<331 
(34) 


The quantities (pM and nP 1 '' can he solved from (31 )*(34) and (11a). with the result as in 
the following expressions: 


„W~ J~, ; t_.. 

2k 1 4 y 3 «+* cofi&o 

(35) 

&TE J 2 4- \/l 4 Cos$ o 

2k 1 A pi ^ cos&o ' 

(3(5) 

It is seen that' the best location for Die monofile is on the negative 
dependent on the angle 4>o of the incoming wave. 

x axis at a distance 
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DIFFRACTION FIELDS 


To see tin effect of the shift of the motiofih* from tin* origin to the complex point, Id 
us consider the diffracted far field and the far fields arising from the monofile located both 
at the edge and in complex space, The far field arising from the exact diffraction current 
(20) can be evaluated in the form 


i/J'M 


, < )k >’ 2 v/tt 

to, ^Mwp’Vj 


n ( ainfy \ 
} \sm|’ + costy ) 


The line current iP 1 at the edge gives rise to the following field: 


eP’ 


p JL -f ( * in 3„ \ 

*' y/tijirtp •/] yJ. 4 cosfy J ' 


(37) 


(38) 


which is seen to coincide with (37) at 4> * tt/ 2, The line current at the point x =* o' ^ 
defined by (35), radiates the field 


1?TM 


E, 


1'. C ' }kP 2y/7T / \ / 1 CQ3<ft \ 

' y/EJWp sfi y ^ + costy J exp \ 2 y /2 J, + costy J ' 


(39) 


By expanding both (37) and (39) at <f> ~ $ 4 6, in power series in terms of f, which 
is a small quantity* it can be seen that the first term i$ coincident with (38), and the 
second terms in (37) and (38) are equal. This shows that (39) approximates (37) more 
closely near the radiation direction <f> as 90°, Examples of direction patterns for these fields 
can be seen in figures 2, 3 for the incidence angles =.60°, 120°. It is seen that while 
the current at the edge only gives omnidirectional radiation, the complex space monofile 
gives quite good estimate for the diffraction patterns for angles about 45° 180°, wheras 

the diffraction field in the direction of the half plane is poorly approximated. From the 
previous we conclude, that the complex space line current gives the diffraction pattern and 
its derivative correctly for the angle 90°. It is obvious that this is the best that can be 
done with an approximation by a single monofile, The value of cou | ( j chosen also 
to meet some other requirement, For example* we might have wanted the correct field in 
the direction 0 ~ 0, for which the coefficient of coho within the exponent function of (39) 
could have been accordingly chosen, In this case, however, it is obvious that the radiation 
in the direction 0 *?■ 180° would be greatly deteriorated. 

Corresponding expressions for the TE field can be written, respectively, 



\gain. alter '*• »i 1 1 « * algebra. if ran In* »•'#*!» that tli* two term approximation* «»f ( tfh him! 
t IJi omim id* lot ^ • A lor Mimll ami the lit**! term ecpial* that of ( II ), I igure* I. 

r * d* ii»»»iiHt i at* that tli* complex spac* motiofile <a«'mdly a current strip with transverM 
current thiwi gives «|tnt #• nood appf oxiniat ion . but only for angles hO".. I Hi r . I hi* i 
1 1 tat ti l\ b* i miim* a flint*- dipolar current cannot pro*lu « a non null radiation field hi its own 
dir* < f ion. w h* n a* t lit* infinite original < m n ut « an. ’I he current strip at tin* edge gives vi *w 
a sinusoidal transverse pattnn. and tli*- complex location distorts that pattern to uiftk* 
th* pattern and it- derivative matrh tha* of th* exact solution 


Ol \l l< M IN( 11)1 N( h 

I lit approximation < an he generalized tot he plane wav** in* ident obliquely tot he »*ug** 

I Ik- v^tvc with I I and I \1 polari/at ions are still considered separately If th** angle with 
th** edge i v the previous theory applies with the dependence i/;ii jkzto.%H n) added 
to th** » mi* lit and fields and with k everywhere in the formulas replaced by ksin^). ft 
p 1-7 Mm t| » normalized line current* I II), (12) and th** location* of the moiiofil* - 
i 'Mi. ( 'll* » must he divided I *\ 

I he incoming plan* wave must be decomposed in two plrrne waves of II and I M 
polari/ation I Ins ran he done with the following expansion for the unit dyadic 

k-l (k > u 2 )[k > u : ) (k (ft ' u : ))(k * (k * ii.) 

P (1 * I 

H*r» . t Ik* first term does not affect w hen multi pi ying • lie fields. miic e k / 0 and k // ft 

Multiplying ( 1'i) by / give us th* decomposition I 0 / ^ * / ‘ Multiplying by 

// w»- hav** tie de* omp*isition // (t • // II 1 ' line, any plane wav* not incident 

along th* edge * an h* decomposed and the monofilar c urrents for each partial field ren he 

summed . 


\< KNOW LKIXiEMKN T 

I hi** vv*#rk was don** while tin- author was on sabbatical leave at the Mil Re^arel 
I.ahoratory of I ilr< tronh *. during th** academic year 19H(. s7, The author f« <|s indebted 
to prof* *s «r .1 . \ Kott* 1* »i this possibility ami fiuitful c!i -»« u» ion- * >on the- pro* i * t-.pic 





Mi I I Nl Mis 


1 \ Smihiim H* M. MmI li# ni«ihs< lir I liroi i*' ilri Dillrai timt. Mailt \tin . \<»l 1 7 . |»|> 

:il 7 it i. i*«m. 

2 I I Knott, l li \ Senior. '*( «»l tin limit Ireiptenev diflra« tion teeli 

iiiijiii v * li* m III I vol. t»2. im* 1 1 p|* I 1*71. \i»\ IM7 I 

.1 I \ LiimIi’II. K I \iki»*kiiM ii. ( *»iii|»|rx N|»a« i mult i|n»lr 1 li»or\ (or k< atterinii prol* 
lfMifc" ’ . I /fs/ lii l rrim IhhihI Siin/» <*n tlnitninaniivtH I In tin, \ kadeiniai hiado. Hu 
dape»t. ltlHti. pp. M I . 

I 1 \. Liiidrll. "Complex spare multipoP- "xpatiHon theory with application to ' fitter 
urn from dielretrii bodies . // / I Iran *> \nirnmt * Pro/ia# . vol AP TV pp. XXX XXX. 
June In 7. 

f » .1 van Bladel." I l»» imilt ipole expansion revisited’ \n h I I* kit I brrtr.' vol. *11. 
pp. 107 III, 1**77. 

ti (i H. Inline (•town run! Ilimrx of Diffraction (oi I. Ire t roniBgiirti* Hair London, 
I upland Peter Peregrinin*, l'*7ti 







PIG CRE CAPTIONS 


Figure I Basic geomet ry of plane wave incidence on a conducting half plane. The angle 
of incidence is bet ween 0 and 180° and the angle $q between 0 and 90 r \ 

Figure 2 TM polarized diffraction field patterns for 0(j ~ 0 and <f ) o « 60°, Solid line: 
exact diffraction field, equation (37), dot line: monofilar current at the edge* equation (38), 
dashed line; monofilar current in complex space, equation (39). The quantify displayed is 
actually the bracketed term in the field expressions. 

Figure 3 Same as Fig.2, but with «/>o = 120°. 

Figure 4 TE polarized diffraction field patterns for == 0 and ^>q ** 90°, Solid line: 
exact diffraction field, equation (40), dot line: monofilar current at the edge, equation (41), 
dashed line: monofilar current in complex space, equation (42). The quantity displayed is 
actually the bracketed term in the field expressions. 

Figure 5 Same as Fig.4, but with <j>Q = 120°. 
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